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Abstract. It is known that the continuous wavelet transform of a function / 
decays very rapidly near the points where / is smooth, while it decays slowly 
near the irregular points. This property allows one to precisely identify the 
singular support of /. However, the continuous wavelet transform is unable to 
provide additional information about the geometry of the singular points. In 
this paper, we introduce a new transform for functions and distributions on R 2 , 
called the Continuous Shearlet Transform. This is defined by SHf(a,s,t) = 
(f,ipast), where the analyzing elements tpast are dilated and translated copies 
of a single generating function if) and, thus, they form an affine system. The 
resulting continuous shear lets ip aB t are smooth functions at continuous scales 
a > 0, locations t G R 2 and oriented along lines of slope s £ R in the frequency 
domain. The Continuous Shearlet Transform transform is able to identify 
not only the location of the singular points of a distribution /, but also the 
orientation of their distributed singularities. As a result, we can use this 
transform to exactly characterize the wavefront set of /. 



It is well-known that if ip is a 'nice' continuous wavelet on R™, and / is a function 
that is smooth apart from a discontinuity at xo S R™, then the continuous wavelet 
transform 



decays rapidly as a — > 0, unless t is near xq I21| . As a consequence, the contin- 
uous wavelet transform is able to resolve the singular support of the distribution /, 
that is, to identify the set of points where / is not regular. However, the transform 
W/(a,i) is unable to provide additional information about the geometry of the set 
of singularities of a more general distribution /. For example, in many situations 
one would like to identify not only the location of a certain distributed singularity, 
but also its orientation. This is very relevant, in particular, in the study of the 
propagation of singularities associated with partial differential equations |17l I23| . 

In this paper, we introduce a new two-dimensional continuous transform, called 
continuous shearlet transform, that is mapping a tempered distribution / 6 <S'(R 2 ) 
to SHf(a,s,t), where a > 0, s £ R and f £ I 2 , are the scale, shear and loca- 
tion variables, respectively. The transform is defined by SHf(a,s,t) — (f,ip a st), 
where the analyzing elements ipast, called continuous shearlets, are dilated and 
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translated copies of a single generating function ip, thereby constituting a contin- 
uous 2-dimensional wavelet-like system. The generator t/j £ L 2 (R 2 ) is chosen to 
be arbitrarily smooth and has compact support in the frequency domain with a 
needlelike structure to capture directionality. For the dilation matrix we employ 
the product of a parabolic scaling matrix associated with a > and a shear matrix 
associated with s £ M., whereas ( e I 2 serves as the translation parameter. The 
system {ip as t : a > 0, s £ R, t £ M. 2 } will be shown to satisfy a Calderon-type 
formula with respect to a special measure, that is, it forms a reproducing system 
ior L 2 (R 2 ). 

We will prove that the continuous shearlet transform is able to identify not only 
the singular support of the distribution /, but also the orientation of distributed 
singularities along curves. In particular, the decay properties of the continuous 
shearlet transform as a — > precisely characterize the wavefront set of / (see 
Section [5J ■ 

Historically, the idea of using wavelet-like transforms to identify the set of singu- 
larities of a distribution can be traced back to the wave packet transform, introduced 
independently by Bros and Iagolnitzer pQ and Cordoba and Fefferman |10) . More 
recently, Smith [221 and Candes and Donoho EI have introduced a continuous 
transform that uses parabolic scaling in polar coordinates and has the ability to 
identify the wavefront set of a distribution. In particular, the Continuous Curvelet 
Transform (CCT) of Candes and Donoho has properties similar to the continuous 
shearlet transform presented in this paper. Namely, in both cases, the wavefront set 
is exactly characterized by the decay properties of the continuous transform. Un- 
like the CCT, however, our approach exploits the framework of affine systems. As 
a consequence, the continuous shearlet transform is not only much closer in spirit 
to the traditional continuous wavelet transform, but also avoids the more compli- 
cated structure of the curvelet construction, which uses infinitely many different 
generators. 

Another motivation for this investigation and the use of the framework of affine 
systems comes from the study of discrete wavelets, and, more specifically their abil- 
ity to approximate efficiently smooth functions with singularities. This property is 
closely related to the micro-localization property of the continuous wavelet trans- 
form. To illustrate this point, let / be a one-dimensional function that is smooth 
apart from a discontinuity at xq and consider the wavelet representation of /: 

where ipj,k(x) = 2~ J/ ' 2 ip(2~°x — k) and ip is a 'nice' wavelet. Since the wavelet 
transform Wf{2~ J ,k) = (f,ipj,k) decays rapidly for j — > oo unless k is near xo, 
it follows that very few coefficients of the wavelet representation are sufficient to 
approximate / accurately. Indeed, the wavelet representation is optimally sparse 
for this type of functions (cf. |2()l Ch.9]). Therefore, the properties of the continuous 
wavelet transform are crucial in order to achieve a sparse representation. 

Traditional wavelets, however, have a very limited capability in dealing with 
discontinuities in higher dimensions. Consider, for example, the wavelet represen- 
tation of a 2-dimensional function that is smooth away from a discontinuity along 
a curve. Because the discontinuity is spatially distributed, it interacts extensively 
with the elements of the wavelet basis, and, thus, many wavelet coefficients are 
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needed to represent / accurately. Related to this is the fact that the corresponding 
two-dimensional continuous wavelet transform is unable to 'track' the discontinuous 
curve. As pointed out by several authors (cf. 01 E]), in order to represent multidi- 
mensional functions efficiently, one has to use representations that are much more 
flexible than traditional wavelets, and have the ability to capture the geometry of 
multidimensional phenomena. One goal of this paper is to show that the study of 
the continuous wavelet transform associated with the affine group on R 2 provides 
a unifying perspective on the construction of such representations. 

The study of the discrete analog of the continuous shearlet transform is cur- 
rently being developed by the authors and their collaborators |121 1131 1141 1151 IT5| . 
In particular, it was shown that, thanks to the mathematical structure of affine 
systems, discrete shearlets are associated with a Multiresolution analysis similar to 
traditional wavelets, which is very relevant for the development of fast algorithmic 
implementations. In addition, shearlets provide optimally sparse representations 
for functions in certain classes and can be easily generalized to higher dimensions. 

The paper is organized as follows. In Section [21 we recall the basic properties of 
affine systems on R™ and the continuous wavelet transform, before introducing the 
continuous shearlet transform in Section|21 In Scction|31we apply this new transform 
to several examples of distributions containing different types of singularities. The 
main result of this paper is proved in Sectional where we show that the continuous 
shearlet transform exactly characterizes the wavefront set of a distribution. Finally, 
in Sectional we discuss several variants and generalizations of our construction. 

2. Affine systems and wavelets 

2.1. One-dimensional Continuous Wavelet Transform. Let Ai be the affine 
group associated with R, consisting of all pairs (a,t), a,tel,(i>0 1 with group 
operation (a,t) ■ (a',i') = (aa',i + at'). The (continuous) affine systems generated 
by ip G L 2 (R) are obtained from the action of the quasi-regular representation ^( a ,i) 
of Ai on L 2 (R), that is 

{4>a,t(x) = 7T( ,t) ip(x) = T t D a ip(x) : (a,t) G Ai}, 

where the translation operator T t is defined by T t ip(x) = %p(x — t) and the dilation 
operator D a is defined by D a ip(x) = a~ 1 / 2 ?p(a~ 1 x). 

It was observed by Calderon 2 that, if ip satisfies the admissibility condition 

(2.1) / |t/>«)| 2 — = 1 fora.e. £ G R, 

Jo a 

then any / G L 2 (R) can be recovered via the reproducing formula: 

/ = / (f,tpa,t)^a,tdfj,(a,t), 

J At 

where d/i(a, t) — dt ^ is the left Haar measure of Ai. Here the Fourier transform 
is defined by ?/;(£) = J ip(x) e~ 2m ^ x dx. As usual, tp will denote the inverse Fourier 
transform. The function tjj is called a continuous wavelet, if ip satisfies (|2.1|l . and 
W/(a,i) = {f,ip a .t) is the continuous wavelet transform of /. We refer to for 
more details about this. 

Discrete affine systems and wavelets are obtained by 'discretizing' appropriately 
the corresponding continuous systems. In fact, by replacing (a, t) G Ai with the 
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discrete set (2- 7 , 2- 7 m), j, m £ Z, one obtains the discrete dyadic affine system 

(2.2) {^ m (x) =T 2im DiiP(x) = DiT m iP(x) : j,m£Z}, 

and ip is called a wavelet if 1)2.2(1 is an orthonormal basis or, more generally, a 
Parseval frame for L 2 (M). 

Recall that a countable collection {-0^}^/ in a Hilbert space 7i is a Parseval frame 
(sometimes called a izgTii frame) for if J^iei l(/> '/'i)! 2 = ll/ll 2 f° r au / S W. This 
is equivalent to the reproducing formula / = Ylieiif^i)^ f° r au / *= ^> where 
the series converges in the norm of 7i. Thus Parseval frames provide basis-like 
representations even though a Parseval frame need not be a basis in general. We 
refer the reader to [HI E] for more details about frames. 

2.2. Higher-dimensional Continuous Wavelet Transform. The natural way 
of extending the theory of affine systems to higher dimensions is by replacing Ai 
with the full affine group of motions on R™, A n , consisting of the pairs (M, t) £ 
GL n (M) x R" with group operation (M, t) ■ (M', if) = {MM 1 , t + Mt'). Similarly to 
the one-dimensional case, the affine systems generated by tp £ L 2 (W l ) are given by 

{fl>M,t(x) = T t D M ^{x) : {M,t) £ A n }, 

where here the dilation operator Dm is defined by Dm ip(x) = | det M\~^ip(M~ 1 x). 
The generalization of the Calderon admissibility condition to higher dimensions and 
the construction of multidimensional wavelets is a far more complex task than the 
corresponding one-dimensional problem, and yet not fully understood. We refer 
to ^3 for more details. 

Now let G be a subset of GL n (R) and define A C A„ by A = {(M,t) : M £ 
G,t£ R n }. If there exists a function ip £ L 2 (R n ) such that, for all / £ L 2 (R n ), we 
have: 

(2.3) f = J^ Jjf,T t D M iP)T t D M iPd\(M)dt, 

where A is a measure on G, then ip is a continuous wavelet with respect to A. The 
following result, that is a simple modification of Theorem 2.1 in gives an exact 
characterization of all those £ L 2 (M. n ) that are continuous wavelets with respect 
to A. The proof of this theorem is reported in the Appendix. 

Theorem 2.1. Equality is valid for all f £ L 2 (W l ) if and only if for all 

CeK"\{0}, 

(2.4) A(V0(f)= / |^(M'0| 2 |detM|dA(Af) = 1. 

JG 

The choice of the measure A on G is not unique. If G is not simply a subset of 
GL n (M), but also a subgroup, then we can use the left Haar measure on G which 
is unique up to a multiplicative constant. Also, observe that Theorem 12.11 extends 
to functions on subspaces of L 2 (M. n ) of the form 

L 2 {VY = {f£ L 2 (R n ) : supp/ C V}. 

2.3. Localization of Wavelets. The decay properties of the functions ipM.t — 
T t Dm ip, where y) £ Cg° , are described by the following proposition. 
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Proposition 2.2. Suppose that ip E L 2 (M. n ) is such that ip € C^°(R), where R — 
suppip C K™. Then, for each k £ N, there is a constant Ck such that, for any 
x S M™, we have 

\i>M,t{x)\ < °k I detM\-i (i + \ M -\x - t)\ 2 y k . 

In particular, Ck — km{R) (||"0||oo + || A fe -0|| oo ), where A = Ym=i ~§s? * s the fre- 
quency domain Laplacian operator and m(R) is the Lebesgue measure of R. 

The proof of this proposition relies on the following known observation, whose 
proof is included for completeness. 

Lemma 2.3. Let g be such that g £ C§°(R), where R c R n is the suppg. Then, 
for each k S N, there is a constant Ck such that for any x S W 1 

\g(x)\<Ck(l + \x\ 2 )- k . 

In particular, Ck = km(R) (||g||oo + ||A fc 5|| 00 ). 

Proof. Since g(x) — j R <?(£) e 27 ™^ d£, then, for every i£l 2 , 

(2.5) \g(x)\ <m(R)\\g || 
An integration by parts shows that 

f Ag(0e 2 ^dt = ~(2n) 2 \x\ 2 g(x) 

JR 

and thus, for every i£l 2 , 

(2.6) (27T\x\) 2k \g(x)\<m(R)\\A k g\\ 00 . 
Using l|2.5f) and l|2.6fl . we have 

(2.7) (1 + (2tt \x\) 2k ) \g(x)\ < m(R) (\\g\U + HA^IU). 
Observe that, for each k 6 N, 

(1 + MY < (1 + (2tt) 2 \x\ 2 ) k < k (1 + (2tt \x\) 2k ). 
Using this last inequality and l|2.7J) . we have that for each x € W 1 
\g{x)\ < km(R) (1 + \x\ 2 )- k (||s|U + HA^U). □ 

A simple re-scaling argument now proves Proposition ^. 21 
Proof of Proposition [2~T2l A direct computation gives: 

il>{M-\x-t)) = [ V^e 2 ™**" 1 ^-')^ 

JR 

= I iPiOe^-^^dt 

JR 

= [ ^{M t n)e 27Tl{x ' t)ri \det M\drj. 

JRM- 1 

It follows that 

\ip(M~ 1 (x -t))\< m{Nr l R) | det M\ ||^(Af*.)||oo = m{R) H^IU 
Using a simple modification of the argument in Lemma 12.31 we have that 
(2tt \M-\x ~ t)\) 2k \ip(M~ 1 (x -t))\< m(R) ||A*^||oo. 
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Next, arguing again as in Lemma 12.31 we have that 

\ip(M~ 1 (x - t))\ < km(R) (1 + \M~\x - t)\ 2 )- k (H^Hoo + IIA^Hoo). 
This completes the proof. □ 

3. Continuous Shearlet Transform 

3.1. Definition. In this paper, we will be interested in the afhne systems obtained 
when A is a subset of A2 of the form 

(3.1) A = {(M,t) : M G G,t £ M 2 }, 

and G c GL 2 (M) is the set of matrices: 

{/a -y^as^ 
M = M as =\ |. a^I.seS 

where / C R + , ScM. It is useful to notice that the matrices M can be factorized 
as M — B A, where B is the shear matrix B — Q "j and A is the diagonal 

/ a \ 

matrix A = .In particular, A produces parabolic scaling, that is, f(Ax) — 

\ y/a J 

X1 W ■ ■ - 2 



X2 



f ( A I ) ) leaves invariant the parabola x\ — x 2 ,. Thus, the matrix M can be 



interpreted as the superposition of parabolic scaling and shear transformation. 

We will now construct a continuous wavelet on R 2 associated with the subset A 
of the affine group, given by (|3.1|) . We will consider two situations, corresponding 
to / = R+ S = R or / = {a : < a < 1}, S = {s G R : |s| < s }, for some s > 0. 

For £ = G i 2 , 6 ^ 0, let V be given by 

(3-3) ^(0=^i.6)=lte)&@- 

Proposition 3.1. Let A &e #wen 6?/ f ff. J|) and fff.jjj) mt/i I = R+, S = R, and 

r/» G L 2 (R 2 ) &e #ttien &y where: 

(i) i/'i € i 2 (R) satisfies the Calderdn condition \2.1\l : 
(n)\m\v = l- 
Then ip is a continuous wavelet for L 2 (R 2 ) with respect to A. 

Proof. A direct computation shows that M*(£i, £2)* = (a£i, a 1 ^ 2 (^2 — s £i))*- By 
choosing as measure d\(M) — , dc f a M ^ ds, the admissibility condition 12. 4|) becomes 

(3.4) A(^)(0=// |^i(aCi)| 2 |^ 2 (a-5(|- S ))| 2 a-idads = l. 

JR JR+ 

Thus, by Theorem 12.11 to show that ip is a continuous wavelet it is sufficient to 
show that 1)3. 4|l is satisfied. Using the assumption on ip\ and ^2, we have: 



AW>)(0 = / / |^iKi)| 2 |<Ma~(t-s))| 2 a-*dacfa 

JR JK+ 

hMa£i)| 2 (/ |^(a^|- S )| 2 d S )^ 



/ l^i(«Ci)| 2 — = 1 fora.e. GR 2 . 
Jr+ a 
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This shows that equality (|3.4|1 is satisfied and, hence, ifi is a continuous wavelet. 
□ 

If the set S is not all of R, then we need some additional assumptions on ip. 
Consider the subspace of L 2 (R 2 ) given by L 2 (C) V = {/ G L 2 (R 2 ) : supp/ C C*}, 
where 

C = {(6,6)eK 2 :|ei|>2and|||<l}. 
We have the following result. 

Proposition 3.2. Let A be given by and with I = {a : < a < 1}, 

S = {s € R : |s| < 2}, and ip G L 2 (R 2 ) be given by fO)) w/iere: 

(i) -01 S i 2 (K) satisfies the Calderbn condition \2. and suppi/ji C [—2, — |]U 

(ii) HV^Hl 2 = 1 suppi>2 c [—1,1]. 

T/ien ip is a continuous wavelet for L 2 (C) V n/ii/i respect to A, that is, for all f G 
£ 2 (C)\ 

= / / / (f,1past)lpast(x)—^dsdt. 
Jr2 7-2 Jo a 

Proof. We apply again Theorem 12.11 for functions on L 2 (C) V . Using the as- 
sumptions on ip2, S and / we have that, for £ G C: 

^(fr+ 2 ) . r 1 . 

|V 2 (s)| 2 rf.s= / |V 2 (s)| 2 ds = l. 

7s(ff- 2 ) ^ 

Thus, for a.e. £ G C we have that 

A(V)(0 - T f \i>{MU)\ 2 a-idads 

J -2 JO 

2 / IV^iKi)! 2 |^2(a-'(|f - S ))| 2 a~i dads 

-2 JO 

1 l ^ £2. _|_2) 

l^iKi)l 2 P ?1 \Ms)\ 2 ds- 
V(a&)| 2 -. 

a 

Since £1 > 2, using the assumptions on the support of ^1 and condition Q2.1f) . from 
the last expression we have that, for a.e. £ G C, 

AW(0 = t X \Ma)\ 2 v = f l<M«)| 2 - = r l^^l 2 - = L 
Jo « J| a Jo a 

This shows that the admissibility condition 12.4fl for this system is satisfied and 
this completes the proof. □ 

There are several examples of functions ipi and ip2 satisfying the assumptions of 
Proposition 13 . II as well as Proposition 13. 21 In addition, we can choose tpi, such 
that ipi, ip2 G (see ^JEl for the construction of these functions). 

Now we can define the continuous shearlet transform: 
Definition 3.3. Let ip G L 2 (R 2 ) be given by lpH)l where: 
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(i) ipi G L 2 (R) satisfies the Calderon condition (f!TT|) . and xpi G Cg°(M) with 
supp^i C [-2,-J]U[§,2]; 

(ii) \\1p2h' = 1, and?/> 2 G Cg°(R) withsupp^ C [-1,1] and i> 2 > on (-1,1). 
The functions generated by under the action of A, namely: 

4>ast(x) =T t D M ^(x) = a"*V n°"j! S j (x-t)j , 

where a G / C R + , seScR and ( £ I 2 , are called continuous shearlets. The 
continuous shearlet transform of / is defined by 

SH f (a,s,t) = (f,ip ast ), ae I cR+,s G 5 Cl,i el 2 . 

Observe that, unlike the traditional wavelet transform which depends only on 
scale and translation, the shearlet transform is a function of three variables, that 
is, the scale a, the shear s and the translation t. Many properties of the continuous 
shearlets are more evident in the frequency domain. A direct computation shows 
that 



Thus, each function ip as t is supported in the set: 

suppvw c {(a, 6) : a e [-§,-£] u [£, f], ig - s\ < 

As illustrated in Figure 1, each continuous shearlet ip as t has frequency support on 
a pair of trapezoids, symmetric with respect to the origin, oriented along a line of 
slope s. The support becomes increasingly thin as a — > 0. 

When 5 = M. and J = R + , by Proposition l3.il the continuous shearlet transform 
provides a reproducing formula Ij2.3[l for all / G L 2 (M. 2 ): 



[ r r \SH f (a,s,t)\ 2 ^dsdt. 
Jr 2 J-oc Jo a 
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On the other hand, if S, I are bounded sets, by Proposition 13.21 the continuous 
shearlet transform provides a reproducing formula only for functions in a proper 
subspace of L 2 (R 2 ). However, even when S, I are bounded, it is possible to obtain 
a reproducing formula for all / G L 2 (R 2 ) as follows. Let 

V3<%) = v3 ( %i,6) = <M6) ^ 2 (|), 

where t/>i, %ji 2 are defined as in Definition liOl and let A (t,) = {(M,t) : M G G {v) ,t G 
R 2 }, where 

(3.5) G {v) = \m = M as = ( ] , a G J, s G s\ , 

[ \-y/as a) J 

Then, proceeding as above, it is easy to show that i})( v ' is a continuous wavelet for 
L 2 {C^"'Y with respect to k^ v \ where is the vertical cone: 

C (,;) ={(6,6)eK 2 :|e 2 |>2and|||>l}. 

Accordingly, we introduce the shearlets ip a v J t = T t Dm^j (v \ for (M,t) G A^ v \ and 
the associated continuous shearlet transform SHf (a, s,t) = (f, Wast)- Finally, let 
W(x) be such that W(£) G C*°°(M 2 ) and 

(3.6) |^(0| 2 + Xc7 1 (e) / 1 |V-iKi)| 2 -+Xc 2 (0 / 1 |^i(a6)| 2 - = l, 

Jo a Jo a 

for a. e. £ G M 2 , where d = G « 2 : ||l < 1}, C 2 = {(6, 6) G 

R 2 : 1 > 1}. Then it follows that W is a C°° window function in M 2 with 
#(0 = 1 for £ G [-1/2, 1/2] 2 , = outside the box {£ G [-2,2] 2 }. Finally, 

let (P Cl f) A = fxci and {P C J) A = f Xc 2 - Then, for each / G L 2 (M 2 ) we have: 

Il/H 2 = [ \(f,T t W)\ 2 dt+ [ f f\SH (Pcif) {a lSl t)\ 2 ^d S dt 

(3.7) + IJ^\^l 2f) (a,s, t )\ 2 ^ dsd t. 

The proof of this equality is reported in the Appendix. Equation (|3.7|l shows that 
/ is continuously reproduced by using isotropic window functions at coarse scales, 
and two sets of continuous shearlets at fine scales: one set corresponding to the 
horizontal cone C (in the frequency domain) and another set corresponding to 
the vertical cone C^ v \ The advantage of this construction, with respect to the 
simpler one where S = K, is that in this case the set S associated with the shear 
variable is the closed interval S = {s : \s\ < 2}. This property will be important in 
Subsection 14.41 and Section 03 

There are other choices of the subset A, given by (|3.1|) . generating affine systems 
with properties similar to the continuous shearlets. Variants and generalizations of 
this construction will be discussed in Sectional 
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3.2. Localization of Shearlets. Since the continuous wavelet ip associated with 
the continuous shear let transform satisfies tf) G C^°(R 2 ), it follows that the contin- 
uous shearlets decay rapidly as \x\ — > oo, that is: 

ipast(x) — 0(\x\~ k ) as |x| — > oo, for every k > 0. 

More precisely, we have the following result. 

Proposition 3.4. Let ip G L 2 (R 2 ) be the continuous wavelet associated with the 
continuous shearlet transform, and let M be defined as in H3.2jl . Then, for each 
k G N, there is a constant C k such that, for any x G R 2 ; we have 

\i> ast (x)\ < C k \detM\-i(l + \M- x (x-t)\ 2 )- k 

= C k a-i (1 + a- 2 {x 1 - h) 2 + 2a- 2 s{ Xl - £i)(x 2 - f 2 ) 
+a- 1 (l + a- 1 S 2 )( 2;2 -i 2 ) 2 )- fc . 

In particular, C'k = k ^ (||^||oo + ||^ fe V'||oo); where A = + is the 
frequency domain Laplacian operator. 

Proof. Observe that, for t = (jj) and x = ( %\ ) in R 2 , we have: 

V>„ flt a = detM-^M" 1 ^-*) = a -4^( V1 _V V2 21 . 

\ a 2 (^x 2 — t 2 j / 

The proof then follows from Proposition ^. 21 where 

R = {(6,6) : 6 G h|, -£] U [£, f], |* - || < v^}. 
It is easy to check that m(R) — ^ ■ D 

4. Analysis of singularities 

As observed above, the mother shearlet ip, constructed in Section satisfies 
ijj G Cg°(R 2 ). It follows that ip G 6>(R 2 ) and, therefore, the continuous shearlet 
transform SHf(a,s,t) — (f,ip as t), a > 0, s G M,< G R 2 , is well defined for all 
tempered distributions / E <S'. 

In the following, we will examine the behavior of the continuous shearlet trans- 
form of several distributions containing different types of singularities. This will 
be useful to illustrate the basic properties of the shearlet transform, before stat- 
ing a more general result in the next section. Indeed, the rate of decay of the 
continuous shearlet transform exactly describes the location and orientation of the 
singularities. 

In order to state our results, it will be useful to introduce the following notation to 
distinguish between the following two different behaviors of the continuous shearlet 
transform. 

Definition 4.1. Let / be a distribution on R 2 , SHf(a, s, t) be defined as in Defi- 
nition rOl and let r G R. Then SHf(a, s, t) decays rapidly as a — > 0, if 

SHf(a,s,t) = 0(a k ) as a — > 0, for every k > 0. 

We use the notation: SHf(a, s, t) ~ a r as a — > 0, if there exist constants < a < 
(3 < oo such that 

aa r < SHf(a, s,t) < (3a r as a — > 0. 
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4.1. Point singularities. We start by examining the decay properties of the con- 
tinuous shearlet transform of the Dirac 5. 

Proposition 4.2. Ift = 0, we have 

SH$(a, s,i) ~ as a — *■ 0. 

In all other cases, SHs(a, s,t) decays rapidly as a — ► 0. 

Proof. For t = we have 

3 3 

(S, ipast) = -0aso(O) = V(0) ~ a~ 3 as a -> 0. 
Next let t ^ 0. Then 

and, by Proposition ^. 41 for each fc G N, we have 

|-0 ast (O)| < C fe a"i(l + aT 2 t 2 + 2aT 2 sM2 + (1 + a" 1 .s 2 )a" 1 ^)" fe . 

Thus, if t 2 7^ 0, then \ip ast (0)\ = 0(a fc - 3 / 4 ) as a -> 0. Otherwise, if t 2 = 0, ti ^ 0, 
then |^orf(0)| = CKa 2 *- 3 / 4 ) as a -> 0. □ 

4.2. Linear singularities. Next we will consider the linear delta distribution 
v p (xx, x-z) = 5(x\ + px 2 ), p € K, defined by 



{v p ,f)= I f(~px 2 ,x 2 )dx 2 . 
Jm 

The following result shows that the continuous shearlet transform precisely deter- 
mines both the position and the orientation of the linear singularity, in the sense 
that the transform SH Vp (a, s,t) always decays rapidly as a — > except when t is on 
the singularity and s — p, i.e., the direction perpendicular to the singularity. 

Proposition 4.3. Ifti = —pt 2 and s = p, we have 

SH Up {a, s,i) ~ as a — > 0. 

In all other cases, SH Up {a, s,t) decays rapidly as a — > 0. 

Proof. The following heuristic argument gives 

£p(£b&) = / I S{xi+px 2 )e- 2m ^ x dx 2 dxi 



That is, the Fourier transform of the linear delta on M 2 is another linear delta on 
R 2 , where the slope — ^ is replaced by the slope p. A direct computation gives: 



3 



/ ^(a,a _ ^Ki-a _ ^a)e- 27rla_1 « l( * 1+f>t2) ^i 
a" 1 / ^i(Ci)V'2(a"^(p-5)) e - 27rM "« 1 ( tl+pt2 )da 
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Figure 2. Decay properties of the continuous shearlet transform SH XV (a, s, t). 

= a~i Vi 2 (a"5(p- s)) ipxi-a,- 1 {h +pt 2 )). 
If s =/= p, then there exists some a > such that |p — s| > y/a. This implies that 
V>2(o _1 ' 2 (p — s)) = 0, and so (z> p , ?/> as t) = 0. On the other hand, if t\ = —pt 2 and 
s=p, then ip 2 (a- 1 / 2 (p~ s)) = $ 2 (0) 7^ 0, and 

ifipAast) = ^2(a~^(p- s)) V»i(0) ~ as a _> 0. 
If ti 7^ — pt 2 , by Proposition 12. 21 we observe that, for all k e N, 

< a-itp 2 {a-^(p-s)) |Vi(a _1 (ii +p*a))| 

< C , fea"^2(a _ 5(p- s ))(l + a - 2 (t 1+p t 2 ) 2 )- fc = 0(a 2 ' c -3) as a -> 0. □ 

4.3. Polygonal singularities. Here we consider the characteristic function \v of 
the cone V — {{x\,x 2 ) : x\ > 0, 5x1 < ^2 < P^i}, where < q < p < 00. We have 
the following result. 

Proposition 4.4. For t = 0, if s = — ^ or s = — i, we have 

SH Xv (a, s, t) ~ a 3 as a — > 0, 
and i/s 7^ and s ^ - we Ziave 

SH Xv (a, s, i) ~ a 1 as a — > 0. 
Fort 7^ 0, if s — —- or s = —- , we Ziaue 

SH Xv (a, s, t) ~ a 3 as a — * 0. 
7n all other cases, SH Xv (a, s, t) decays rapidly as a — > 0. 

The decay of the continuous shearlet transform of xv is illustrated in Figure 2. 
As shown in the figure, the decay oiSH Xv (a, s, i) exactly identifies the location and 
orientation of the singularities. It is interesting to notice that the orientation of the 
linear singularities can even be detected considering only the 'point singularity' at 
the origin. 
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Proof of Proposition HTH The Fourier transform of xv can be computed to be 
XV(£i,6) = Cj— — l ( — —t, where C = ^ + ^ . 

(6. + «6)(6. +K2) (27r) 2 

A direct computation gives: 

= CJ I / 77-TWT7T^Ki)^( a ^(| - S ))e 2 "^6rf6 
Ca"^ / / , — -i — — & (£1)^(0-^0^ 



(a ^i + ?6)( a 6 +K2) u 
. e a»*(«~ 1 Ci.&)*df 1 df 2 

^ 4 11 (a-^i +9£i(«- 1/2 6 + a-i S ))( a -i6 + + a^s)) 

•^i(£i) ^2(6) e 2 " ?l(a " (tl+st2)+a " /2 « 2t2) d£i d£ 2 

^'11 (a- 1 / 2 ^!! + sq) + qji J(^ 1/2 Ci(l + *P) 
•^i(a) ^2(6) e 2 -^(-~ 1 (*i+ s ^)+ a ^ 1/2 € 2 t 2 ) da ^ 

Let us first consider the case t = 0. By the previous computation we can rewrite 

(Xv,i>aso) as 

Ca * I I 1 -i/2c n 1 — \ T TTYf -i/2c n _l — \_l f t ^ ife)^fe)^4 
If s 7^ — i and s 7^ — ~, for a<C 1 we can rewrite (xv, 4>aso) as 

Cai I I TFTTZ \~i — 1/2 c £ ut n i r~i — 1/2 £ c \ ^1(^1)^2(6)^1^2 

JrJr (?i(1 + «?) + a ' Ki6)(si(l + «P) + a 7 K1C2) 

~ Cat / / ^ ffi 

•Wr (£i(1 + «<?)) (£i(l + sp)) 

hence 

(xv, V'aso) ~ a 3 as a -> 0. 
The above computation also shows that if t = and s = or s = -i, we have 

(xv, V'aso) ~ a 3 as a — * 0. 

Next, let us consider the situation, where t lies on one singularity but t =/= 0, 
i.e., ti = pti or ti = qt\. Here we will only examine the first case. The second 
one can be treated similarly. First let s = — i, i.e., s is perpendicular to the linear 
boundary of the cone X2 — px\. For a«l we have 

(xvAast) = Cai f f §L ^1(6)^2(6) 

JrJr {a 1/2 Ci(l - Q/P) + <?£i6)Ki6 



= Cat 



(a(i-g/p) + a 1/2 g66Ki6 

•^i(6)^2fe)e 2ma " /2ptl?1?2 ^irf6 

3 

a 4 as a — > 0. 
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Secondly, let s ^ — ^. We have: 

(XV,^ast) 

Si 



Ca~ 



(o-V2^(l + sq) + gei6)(a- 1 / 2 ei(l + sp) + a-Vap&fc) 



= Ca* / ^(eO^i^Oe 2 " " 1 * 1 ^ 8 ^ 1 ^!, 

where 

^ = 1 (FVa£(I + sg) + ^2X^6(1 + s P ) + p£i6)^ (6) 

,g27ria _1/2 tipfi^ 2 ^ 2 

Since ^1 and ^2 are band- limited, the function tp has compact support, hence 
{tpi^iY i s °f rapid decay towards infinity. Thus 

(Xv^ast) = Ca*{<p$i) y {arHi{l - sp)) = 0(a k ) as a -► 0. 

Finally, in case £2 7^ P^ij ^2 7^ 9^1 and t\ 7^ 0, a similar argument to the one above 
shows that {xv,^ast) decays rapidly also in this case. □ 

4.4. Curvilinear singularities. We will now examine the behavior of the contin- 
uous shearlet transform of a distribution having a discontinuity along a curve. 

Let B{x\,X2) — xd{xi,X2), where D — {(acijXa) G ^ 2 : x \ + x \ — !}■ We have 
the following: 

Proposition 4.5. If t\ + t\ = 1 and s = j-, t\ 7^ 0, we /iave 

SHb{o-i s, t) ~ a 3 as(n 0. 
In aH other cases, SHb{o-, s, t) decays rapidly as a — ► 0. 

The assumption 7^ shows that the shearlet transform SHb(u, s,t) is unable 
to handle the vertical direction s — > 00. To provide a complete analysis of the 
singularities of B, we need to use both SHb(cl, s,t) and SH% ''(a,s,t) (as defined 
in Section |3J). Since the shear lets tp^ st are defined on the vertical cone using 
one can obtain a similar result to Proposition ^. 51 for s = I 1 , t% 7^ 0. 
Since the argument for both cases is exactly the same, we will only examine the 
transform SHb(o,, s, t). 

In order to prove Proposition 14.51 we need to recall the following facts. First, 
we recall the asymptotic behavior of Bessel functions, that is given by the following 
lemma (cf. |24j ) : 

Lemma 4.6. There are constants Ci,C% such that 

Ji(2ttA) ~ \~i (Ci e 2mX + C 2 e~ 2mX ) as A 00, 
and, /or N — 1,2, ... , constants C± , C2 such that 
d \ N 



(d\) 



-- ) Ji(2ttA) ~ A" 5 (Cf e 2 " A + Cf e - 2 ™ A ) as A 



Secondly, we recall the following fact concerning oscillatory integrals of the First 
Kind, that can be found in |^| Ch.8]: 
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Lemma 4.7. Let A e C£°(R) and $ G C 1 ^), with $'(*) ^0 on supp A TTien 
J(A) = / A(t)e 2mX ^dt = j;" 1 ^ / r>"(i4(t)) e 2 ^*Wrft, 

il (2mA) J R 

/or iV= 1,2,..., urf»ereZ>(A(i)) = 

We can now prove Proposition 14.51 
Proof of Proposition 14.51 The continuous shearlet transform of B(x) is given 
by: 
(4.1) 

SH B (a,s,t) = {B,% st ) = a i f f V>i«i) - *)) e" 27rie ' B(0 d£i d&. 

JR JR 

The Fourier transform -B(£i,£ 2 ) is the radial function: 

where J\ is the Bessel function of order 1. Therefore, the asymptotic behavior of 
5(A) follows from Lemma f4. 61 with the factor A -1 / 2 replaced by A -3 / 2 . 

Because of the radial symmetry, it is convenient to convert l|4.1|l into polar 
coordinates: 

SHb((i, s, t) 

= ai J J ^ 1 (apcos9) ^ 2 (a~' (tanfl - s)) e" 2 ™^ 1 cos9 +* 2 sine > B(p) pdpd9 

(4.2) =ari J J Mpcos9)^j 2 (a-i{tane - s))e- 27r ^ {tlCOs9+t2Sin9) B(^) pdpdO. 

We will now examine the asymptotic decay of the function SHb{o., s, t) along the 
curve dB for a — ► 0. Thus, we set tf+t 2 = 1 and, without loss of generality, assume 
a < 1. As we will show, the decay will depend on whether the direction associated 
with s is normal to the curve dB or not. 

Let us begin by considering the non-normal case s ^ t^/ti- From (|4.2() . we have: 

SH B (a,s,t) = a-i J I(a, p) pdp, 

where (using the conditions on the support of ipz) 

I(a,p)= [ ^i(pcos6»)^ 2 (a"5(tan6i-s))e- 27rl »( tlCOse+t2Sin ^d6l. 

J | tan 0— s\<,/a 

Observe that the domain of integration is the cone |tan# — s\ < \J~a about the 
direction tan# = s, with a < 1. This implies that 9 ranges over an interval. Since 
the conditions on the support of tpi implies that |pcos#| C 2], it follows that p 
also ranges over an interval and, as a consequence, I(a,p) is compactly supported 
in p. 

We will show that I(a, p) is an oscillatory integral of the First Kind that de- 
cays rapidly for a — > for each p. To show that this is the case, we will apply 
Lemma |4~7I to I(a,p), where A(0; p) = ipx(pcosd) ?/> 2 (a~ 1//2 (tan# - s )), $(9;p) = 
p(ti cos9 + t2 sin 9) and A = a -1 and p is a fixed parameter. Observe that <&'(#; p) — 
p(-tx sin + t 2 cos 9) and ^ for taxi 6 ^ fa. Thus, we have that, if 
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\s — I 3 - j > y/a, then the function $'(0;p) ^ on suppA A direct computation 
gives 



D(A(6-pj) - 9 ^(P cos0 )^( a 2 (tan6»-s)) 



89 p(—tism0 + t2cos6) 

«in fi ^ i 1 

-0i (pcosO) ip2 (a~ 3 (tan 9 - a)) + 



ii sintf — *2 cos 
_i sec 2 



p(t2 cos 6 — ti sin 
i 2 sin + ti cos 



-■01 (p cos 0) V4(a 2 (tan#-s)) 
7/'i(pcos0)V'2(a _5 (tan0 — s)). 



p 2 (ia cos 9 — 1 1 sin 9) 2 

Thus, using the assumptions on ^i> ^2) we have: 

\D(A(0;p)) \ < a"* C(9,p) (||&&l|oo + + Halloo)- 

As observed above, the assumptions on the support of tpi, ip2 imply that D[A{9\ p)) 
is compactly supported in p away from p = 0. Using this observation and $'(0) 7^ 0, 
we have that 

||£>(^)lloc < Ca~i. 
Applying the same estimate repeatedly, we have 

||I?*(A)|U<CWa-3. 
Thus, using Lemma 14. 71 with A = a -1 , we conclude that 

N 

sup \I(a, p)\ < C a" 5 " . 
p 

This implies that, under the assumption that we made for t — (iijia) and s, the 
function SHb{o., s, t) decays rapidly for a — ► 0. 

Let us now consider the function |(.B, where t\ + i 2 . = 1 and s = ^Ai 

(corresponding to the direction normal to 8B). For simplicity, let (ti,t%) = (1,0). 
The general case follows using a similar argument. From 1)4.2(1 . using the change of 
variables it = a -1 / 2 sin0, we obtain 

(4-3) (B, ^ a0(m ) = a"! / B(f ) ?7 a (p) e 2 ^pdp 

where 

Va( P )= [ il+a) ' MpV^^)M t= = J e 2 "^-" ''" 



^-(l+a)- 1 / 2 VT — a« 2 vl — aw 2 

The assumptions on the support of 02 imply that |~/f=^f I < 1- This is equivalent 
to \u\ < (1 + a) 1 / 2 . Similarly, the assumptions on the support of 0i imply that 
\py/l — au 2 \ C [5,2] and, thus, p ranges over a closed interval. As a consequence, 
the functions rj a (p) are compactly supported. For < a < 1, the functions 

e 2jri£(Vl-au 2 -l) 



/l a (u) = ^i(p\/l - au 2 )^ 2 {—E „ ) ■ 



Vl — au 2 \/l — au 2 
are equicontinuous and they converge uniformly: 

lim h a {u) = h (u) = Mp) Mu) e -^ u2 . 
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Thus, we have the uniform limit: 

Urn Va (p) = Vo (p) = [ & (p) i, 2 (u) e-^ 2 dp, 
J-i 

and the same convergence holds for all ^-derivatives. In particular, ||?7a||oo < C, for 
all a < 1. 

Using the asymptotic estimate given by Lemma 14 . 61 into (|4.3|) . we have: 

KA4o(i,o))| ~ a" f (Ci J^Vai^e 4 ^ pdp + C 2 J r, a (p) p dp) 

= ai (Ci F a (-^j +C 2 J F a (p) dp) , 

where F a (p) — T) a (p) p" 1 / 2 . The family of functions {F a : < a < 1} has all its p 
derivatives bounded uniformly in a, and so F a (—^) decays rapidly as a — > 0. On 
the other hand, J F a (p) dp tends to J f]o(p) p^ 1 ^ 2 dp, and thus, combining the two 
terms, we obtain that 

1(^,^0(1,0)) I ~ oJ asa^O. 

If t is not on dB, then one can show that SHb(cl, s, t) has fast decay. This follows 
from the general analysis given in Section [5] □ 

5. Characterization of the wavefront set using the shearlet 

transform 

The examples described in Section 0] suggest that the set of singularities of a 
distribution on M 2 can be characterized using the continuous shearlet transform. 
In this section, we will show that this is indeed the case. In order to do this, it will 
be useful to introduce the notions of singular support and wavefront set. 

For a distribution u, we say that x £ R 2 is a regular point of u if there is G 
Co°(U x ), where U x is a neighborhood of x and 4>{x) ^ 0, such that <pu £ C^°(R"). 
Recall that the condition 4>u € is equivalent to (<pu) A being rapidly decreasing. 
The complement of the regular points of u is called the singular support of u and is 
denoted by sing supp (u) . It is easy to see that the singular support of u is a closed 
subset of supp (u). 

The wavefront set of u consists of certain (x, A) € K 2 x M, with x E sing supp (it). 
For a distribution u, a point (x, A) £ R 2 x K is a regular directed point for u if there 
are neighborhoods U x of x and V\ of A, and a function £ C^°(M 2 ), with cj> = 1 on 
U x , so that, for each N > there is a constant Cjq with 

\(u4>nv)\<c N (i+\r,\r N , 

for all r\ — (771, 772) € R 2 with ^ 6 V\- The complement in R 2 x R of the regular 
directed points for u is called the wavefront set of u and is denoted by WF(u). 
Thus, the singular support is measuring the location of the singularities and A is 
measuring the direction perpendicular to the singularity. 1 

In the examples presented in Section one can verify the following: 

^This definition is consistent with where the direction of the singularity is described by the 
angle 9. Observe that our approach does not distinguish between 9 and 9 + it, since the continuous 
shearlets have frequency support that is symmetric with respect to the origin. However, in Section 
l6lwe discuss a variant of the continuous shearlet transform, which can distinguish these cases. 
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(i) Point Singularity S(x): sing supp (<5) = {0} and WF(S) — {0} x M. 

(ii) Linear Singularity v p (x): singsupp (v p ) = {(— px2,X2) ■ X2 G K} and 
WF(v p ) — {((—px 2 ,X2),p) : x 2 G K}. 

(iii) Curvilinear Singularity B(x): singsupp(-B) = {(11,2:2) : x 2 + x^ — 1} and 
^(5) - {((*!, a*), A) : x\ + x\ = 1, A = a }. 

As observed in Section all these sets are exactly identified by the decay prop- 
erties of the continuous shearlet transform. Indeed, we have the following general 
result: 

Theorem 5.1. (i) Let K = {t Q G M 2 : for t in a neighborhood U of to, 
\SH f (a,s,t)\ = 0{a k ) and \SH { f {a,s,t)\ = 0{a k ) as a^O, for all fc G N, 
with the 0(-) -terms uniform over (s,t) G [—1,1] x U}. Then 

sing supp(f) c = 1Z. 

(ii) Let T> = V x U V 2 , where T>\ = {(t ,So) G R 2 x [-1,1] : for (s,t) in a 
neighborhood U of (so,*o), \SHf(a, s, t)| = 0(a) as a — + 0, for all fc G N, 
with the 0(-)-term uniform over (s,t) G U} and T>2 = {(to, So) G M 2 x 
[1, oo ) : for (i, t) in a neighborhood U of (sq, to), \SHj^ (a, s, t)| = 0(a k ) as 
a — > 0, /or aZ/ fc G N, luit/i t/ie 0(-)-term uniform over (-, t) G Z7}. TTien 

W(/) c = P. 

The statement (ii) of the theorem shows that the continuous shearlet transform 
SHf(a,s,t) identifies the wavefront set for directions s such that |s| = |g| < 1 (in 

the frequency domain). The continuous shearlet transform SH^\a, s,t) identifies 
the wavefront set for directions s such that |s| = 1 1^ | < 1, corresponding to 1 |^-| > 1. 
The proof of Theorem 15 . II will require several lemmata and will adapt several ideas 
from The following lemma shows that if t is outside the support of a function 
<?, then the continuous shearlet transform of g decays rapidly as a — > 0. 

Lemma 5.2. Let g G L 2 (R 2 ) with ||.g||oo < oo. If supp(g) c8c R 2 , then for all 
k>l, 

\SH g (a, s,t)\ = |( 5 , i> ast )\ < C k C(s) 2 || 5 |U a 1 * (l + C(s)- 1 a - 1 d(t, B) 2 y k , 

where C(s) = ^1 + 4^- + (^s 2 + ~\ and Ck is as in Proposition \3.4\ 

Proof. Since |jg||oo < oo, by Proposition 13.41 for all k G N, there is a Ck > 
such that: 

ast)\ S: ||.£?l|oo 

JB 

< Ck\\g\\ooa-i { (l+\\M-\x-t)\\ 2 ) h dx 



(5.1) = C fc HflHoo a-5 f (l+HM-^ll 2 ) 

■>B+t v 



-fc 

dx, 



where M = ^~^ S J. Observe that ||x|| = IIMM" 1 ^ < || M || op \M~ x x\ , and, 
thus, 

\\M^x\\> ' " 



M\\ op 
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Since M = (J (o and || (« Hop = V*, then 

IIM-^H > (7(a)- 1 o - * ||x||, 
where C(s) = || (J I s ) Hop- Using these observation in (|5.1|l . we have that 

|(fl,^ M t)| < c fc Halloo a-5 / (i + c( s )- 2 a- 1 \\x\\ 2 Y k dx 

JB+t K ' 
/>oo 

= C k C{sf\\g\Uofl \ (l + r 2 )- k rdr 

JC(s)- 1 a- 1 / 2 d(t,B) 

= C k C(a) 2 \\gWooai (l + CW-V 1 ^^) 2 )"'. 

At last, we compute (7(s). We have (J ~ s ) (i s °) = Pif "/V The largest 

i 

eigenvalue of this matrix is A max = 1 + %- + ^s 2 + ^-J . Thus we have C(s) = 

(l + ^ + (s 2 + ^) 2 V for all sei. □ 

We can now prove the following inclusions. 

Proposition 5.3. LetlZ andT> be defined as in Theorem \5.1\ Then: 

(i) sing supp{f) c C 1Z. 

(ii) WF(f) c C D. 

Proof, (i) Let to be a regular point of /. Then there exists cf) G C Q yo (M 2 ) with 
<fr(to) = 1 on B(to,S), which is the ball centered at to with radius 5, such that 
(^/ G C°°(M 2 ). We will show that i G 72.. For this, we decompose SHf(a, s, t) as 

(5.2) SH f (a, s, t) = (i> aat ,<j>f) + Wast, (1 ~ 4>)f). 

Observe that we have 

\^ast,<f>f)\ < « f J ^ l^l«i)| \M^( ( i ~ *))| 

In the following, we will estimate the above integral for £i > 0. The case £i < is 
similar. 

Since <j> G C£°(R 2 ), for each k G N there exists a constant Cfe with \4>f{i)\ < 
Cfc |£| _2fc - Using this fact, together the assumptions on the support of f/> os t, for 
k > 2, the first term on the RHS of (|5.2(1 can be estimated as follows: 

at / |&(o&)||&(-W£-«))ll*/(OI# 

< CfeiHiooa* / / ier 2fe rf6^i 

/■! /-(s+aA)6 

< c^iiviiooo^ / / e 2 _fc «2^i 

= Cfc2 ^°° af ((- + v^) 1 - - (- - /* r 2fe dex 
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Cfc 2 ~ fc 1 1 -0 1 1 oo a 4 _ ^i_fc^2-2fc 
S Jfe(2Jfe-l) lV j U ' 

Thus, the above quantity is a 0(a k ) as a — > 0, uniformly over (t, s) g -B(to, I) x K. 

Using Lemma [5.21 we have the following estimate for the second term on the 
RHS of (l5~2l) : 

KV„ rt> (1 - <f>)f)\ < C k C{s) 2 ||(1 - 0)/|U «i (1 + C( S )~ 1 a- 1 d(t, i?(io,^) c ) 2 )- fc , 
where fc g N is arbitrary. Since |j(l — 0)/|j oo <oo and s is bounded, this yields 

uniformly over (i, s) g B(to, |) x [— 1, 1]. A similar estimate hold when SHf(a, s, t) 
is replaced by SHf(a,s,t). This proves (i). 

(ii) Let (to, so) be a regular directed point of /, with s € [— 1, !]• Then there 
exists a <E Cfi°(R 2 ) with </>(io) = 1 on a ball B(to,Si) such that, for each k £ N 
we have | I = + |£|)- fc ) for all £ £ I 2 satisfying f g B(s ,<5 2 ). We will 
prove that (to, s ) g 2?. For this, we decompose SHj(a, s, t) as in l|5.2|) . The second 
term on the RHS of l|5.2[) can be estimated as in the case (i). For the first term 
of H5.2f> . we only need to show that suppVw C {£ £ M 2 : f g -B(s ,<5 2 )} for all 

(s,t) € B(so,£ 2 ) x B(to,6i), since in this cone 4>f decays rapidly. As above, we 
only consider the case £i > 0; the case £i < is similar. The support of ifj ast in 
this half plane is given by 

Let (s,i) g -B(so,<5 2 ) x B(to,Si). The cone {£ g M 2 : f g B{s ,S 2 )} is bounded by 

the lines £ 2 = (so - <5 2 )£i and £ 2 = (so + £2)6.. Now l ct (6,6) € supp^ ast . Then, 
for a small enough, we have 

If -so| < Va~<6 2 , 

and this completes the proof in this case. In case |so| > 1 (this corresponds to 
if I < 1); we proceed exactly as above, using the transform SH^ (a,s,t) rather 
than SHf(a, s, t). □ 

For the converse inclusions we need some additional lemmata. For simplicity of 
notation, in the following proofs the symbols C and Ck are generic constants and 
may vary from expression to expression (in the case of Ck, the constant depends 
on k). 

Lemma 5.4. Let S C K be a compact set, and g g L 2 (M) with ||g||oo < 00. Suppose 
that suppg C B for some B C K 2 and define (B v ) c = {x g M 2 : d(x,B) > 77}. 
Further define h g L 2 (R) by 

POO p P 

h(&= I I I SH g {a,s,t)4> ast (Z)dsdt-^. 



Then /i(£) decays rapidly as |£| — * 00 with constants dependent only on \\g\\oo o,nd 
V- 

Proof. Using the fact that S is compact, Lemma 15.21 implies that, for each 
k > 0, 

\SH g (a, s,t)\< Ck a* (1 + a^d^, B) 2 y k , 
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where Ck depends on ||<?||oo but not on s. By definition, the support of ip a st is 
contained in the set 

(5.3) I>, s) = U G K 2 : | < o|e| < 2, \s - || < 
Thus, |$„<rt(£)l < C"a^Xr(a, s )(0 and 

(5.4) / XT(a,-)(0d8< / r ds<C"^- 

Collecting the above arguments, 

< r [ [ \SH g (a,s,t)\\i> ast (0\dsdt^ 

Jo J(Bn)cJs 01 

< C fc r f /"ax r(O)S) (0(l + a-M*,e) 2 )" fe ^*^ 

JO J{B^)"JS a 

< C k / / X r(a, s )(0^ a (l + a- 1 d(i,e) 2 )- fc ^- 

2 

< C k [™ aT% [ (1 + a- 1 d(i,£) 2 )~ fe dtda 



< C k I aT? / (l + a^r^rdrda 



'i 



1 C k I I'' a~i (l + a-'^y^da 

2TeT 

< c fc ier^(i + ieh 2 r fc+2 - 



Since this holds for each k > 0, /i(£) decays rapidly as |£| — > oo. □ 

Lemma 5.5. Let S" C K and B C K 2 compact sets. Suppose that G(a,s,t) 
decays rapidly as a — > uniformly for (s, t) E S x B. Define h G L 2 (R) &y 

^(0= / / / G(a,s,t)$ ast (Odsdt^. 
Jo Jb Js a 

Then h(^) decays rapidly as |£| — ► oo. 

Proof. As in Lemma [5.41 we will use the fact that |^ a si(£)l < C a ^Xr(a,s)(£)> 
where T(a, s) is given by (|5.3|) . and the estimate 1)5. 4|) . Also, by hypothesis, for each 
k > and a > we have 

sup{|G(o,M)| : If I G [£, £], teB}<C k a k . 
Using all these observation, we have that, for each k > 0, 

\H0\ < r I I \G(a,s,t)\\i> ast (0\dsdt^ 
Jo Jb Js a 

< C k / / Xr(a )S )(£) afc ~ S dsdtda 

Jo JbJs 

< C k i a k ~i da 
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< C k \Z\- k+ i. □ 

The proof of the following lemma adapts several ideas from [3| Lemma 2.3]. 

Lemma 5.6. Suppose < ao < a\ < 1 and |s| < sq. Then for K > I, there is a 
constant Ck, dependent on K only, such that: 

KM,.;)\<c K ( 1 + !lY' ( 1 + \L^.)-' (, . IK'"'')" 



a o / \ a i / \ a i 



Proof. By the properties of tp, we have that, for ||£|| > \ and any k > 0, there 
is a corresponding constant Ck such that 



(i + icii + i^D* 



We also have that = for ||£|| < \. Thus, observing that M* s £ = (a£i, v^6" 
y^s^i), it follows that: 



3 

|^.«t(OI < c * 



(l + a\^\ + V^\^-s^\) k - 

Using polar coordinates, by writing £i = r cos 9 and £2 = r sm8, this expression 
can be written as 

3 

Wast(r, 0)\ < C k — : — r , . Q ~n\Xk~ ' 

(1 + ar\ cos&\ + V ar l sm# — s cos9\j k 
For \9\ < 7r/2, using the assumption |s| < so, the last expression can be controlled 

by 

(5.5) \i> ast (r,e)\ <C k 



k ' 

(1 + a r + ^/ar | sin 6 — s\) 



In addition, since sin# ~ 9 on \9\ < n/2, we can replace sin 8 with 9 in the above 
expression. 

Let As = s — s' , ao = mm(a, a') and ai — max(a, a'). Using (|5.5|) . and applying 
the same argument on |0| < ir/2 and ir/2 < \9\ < it, we have that 

(a a') 3 r 



/ 






/•OC 


<c fe / 










< c 




a = 


sfar 
1+ar' 



l + ar+ x /a7-|6 , -s|) fc (l + a' 
(aa') 3/4 r Z" 00 1 



■ k 

r -\- \fa' r\6 — s 



a- 1 (1 + ar) K ' (1 + a' r) K J-00 (1 + a |0|) fc (1 + a' \9 + As\)' 



d9dr 



■ d8 dr, 



1+0 
1, k > 1, 



/■°° 1 1 

/ 1 p dO < C k T . 

J -00 {1 + a\9\) k {I + a> \9 + As\) k a (1 + a' \ As\) k 
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From the definition a', for r > 1/a', we have that -k -4=? < a' < -4=. Thus, for 
r > 1/a', provided k > 1, the last expression gives 



(5.6) 



-oo (l + a\9\) k (l + a>\9 + As\) k 



:d6 <C k 



(1 +ar) 



\As\ 



Another simple estimate, provided k' > 1, gives: 



(5.7) 



1 



i (l + a rf (l + ai r) k 



dr < C k , - 1 + ^ 



«o 



a i 



«o 



Thus, using (|5.6|) and (|5.7|) . for a > a' > 1, ao = a ' an d a i = a i we have that 



VVs't'(£) Vw(£) 



d£ < c k [ ^ 

a 



a 



1 + — 

a 

-fc+2 



-fc+1 



|A S | 



1 



|A S | 



Similarly, for a > a' > 1, a\= a! and ao = a, a similar calculation gives 

-fc 



VWt'(£) 4>ast(0 



d£ < c k \^ 

a 



a Q 



ao 
-fc+l 



1 - 

|A S | 



|A S | 



In general, renaming the index k, we can show that 



(5.8) 

We have 



4>a's>t>(£) 



d£<C k [l + ^ 
ao 



\As\ 



d - _ d - _ ~ 

^?/w(£) = (tt-Vas)i st (0, -KTlllastiS) = Va^ast{C), 

<9 2 - <9 2 
^TjV'astCC) = (a - \/as) 2 ^ st ((), —2 = a^ast(0- 

Thus, observing that a, a' < 1 and |s| < so, we have that 



A e Vw(OVv s 't'(£) < C'a x |Vw(£)l |Vw*'(£)l- 



Set 



A, 



(2tt) 2 ai 



On the one hand, for each k, we have 



On the other hand 

(5.9) L*( e - 2 **(*-0) = ^ 1 + ' 

Repeated integrations by parts give 

(V>ast,VWt') = / Vw(?)VWt'(£)*; 



-2iri(,(t-t') 
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(a a') 3/4 4(KA) e-^-V d{, 

Therefore, from the last expression, using (|5.8|l - (|5.9|l . we have 

\(^u^)\<c k fi+fc^V* (i + a -±-Y k (i ^ ' As ' N 1 



The proof is completed recalling that, for m > 0, (1 + |x|)~ 2m ~ (1 + |x| 2 )~ m . 
That is, there are constants C x , C 2 > such that d (1 + \x\ 2 ] T m < (1 + \x\)- 2m < 
C 2 {l + \x\ 2 )- m . 

From Lemma 15.61 the following result can be easily deduced. 

Lemma 5.7. Let fa S C°°(R 2 ) supported in £>(0, 1), and define fax) = fa(a7 (x— 
t)). 

(i) Suppose < y^oo < \ZoT < a<p < 1. 27ien /or if > 0, 



(ii) Suppose < y^ao < < ^fal < 1, ai < a^j. Tften /or if > 0, 

|<^, Vw,)l < cx (i + ^) ~* (i + l -^f) " (i + fcf^ 

(hi) Suppose < y/ao < < ai < y'oi < 1. TTiera /or if > 0, 

Now we now complete the proof of Theorem 15. II 
Proof of Theorem 15. U 

Since one direction was proved by Proposition 15.31 we only have to prove the 
inclusions: 

(i) Tl C sing supp(/) c ; 

(ii) V C WF(/) C . 

First we prove (i). Let to € 1Z. Then for all t 6 B(to, 6), a ball centered at to, we 
have that \SHf(a, s,t)\ — 0(a k ) as a — > 0, for all fc € N with the 0(-)-term uniform 

over (t, s) G B(<o, 8) x [—1, 1], A similar estimate holds for SHf"'(a, s, t). 

Choose <fi S C°°(R 2 ) which is supported in a ball -B(io, with i/< { and let 
77 = I . Set g = 0/ and consider the decomposition 

0/(0 = 50 (0 + 51 (0 + 3a(0 + 53(0 +34(0, 

where 5o (0 = (<^(/)) A (0, where P(f) = J R (f,T b W)T b Wdb, W is the window 
function defined by and 

5i(0=XCi(0 / ^a S t(0< SW <?( a jM) rf M a :'M), i = 1, 2, 
•/Qi 
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fli+ 2 (0=Xc 2 (0 / $ ast (QSH g v Ha,s,t)d l i(a,s,t), i=l,2, 

where Ci, C2 are defined after equation dfj.(a,s,t) = ^§ dsdf, Qi = [0, 1] x 

[-2,2] x B(t ,7j) and Q 2 = [0,1] x [-2,2] x B(t ,r]) c . The term <? (£) decays 
rapidly as |£| — » 00 since and P(/) are in C°°. The term 52(C) decays rapidly as 
|£| — > 00 by Lemma l5~4l In addition, by Lemma l5"51 <7i(£) decays rapidly as |£| — > 00 
provided that SH g decays rapidly as a — > uniformly over (i, s) g £?(to, ??) X [—2, 2]. 
We will consider only the analysis of the terms c/i, for i = 1,2; the cases i = 3,4 
are very similar. 

We claim that SH g indeed decays rapidly as a — > uniformly over B(t ,r]) x 
[—2, 2]. In order to prove this, we decompose / as / = P(f) + Pc±f + Pc 2 /, where 
{PcJ) A = fXd and (Pc 2 /) A = /Xcv It is clear that SH$ P (f) decays rapidly by 
the smoothness of </> and P(f). Next, we examine the term Pc ± f- The analysis of 
Pc 2 / is ver y similar and will be omitted. We use the decomposition Pci/ = /1 + /2, 
where 



fi(x) = / ip ast (x)SHf(a,s,t)d(i(a,s,t), £=1,2. 

Let us start by considering the term corresponding to /1. We have: 

■SH^/i (a, s, *) = (<t>fi,ip as t) = / (<pip ast ,ip a > s > t >) SHf(a' , s' ,t') d[i(a' , s' ,t') 

Qi 



We will decompose Qi = Q10 U Qn U Q12, corresponding to a' > <5, a' < <5 < Va 7 
and a/o 7 < <5, respectively. In case i/a, v / a' < ^, by Lemma f5. 71 we have that 



-K , M/ , ,,, ||2 x -K 



(5.10) \{Ma«rfaUv)\<C K fl ' !l</ / ' )| 







(5.11) / + (a')™ 1 777v3 — Cm.if a m , < a < S. 



a J \ ai 
This implies that, for m > 4 and K > to — 1 
r<5 / \ -if 7 / 

a l \ / r\m aa ^ r< - m-2 

a J (a) (a') 
In fact, for a' = ao < a = 01, 

f 1 + 4) ~ K («') m = « m ~ 2 /* a + x )~ K dx = ^ « r ' 

\ a / (a j Jo 

For a = ao < a' = ai, 



1 + - (a') m 7-7^ = a m -" / x m - A (1 + a?) - * dx 



a J x " ' (a') 3 

< a" 1 - 2 I x m -' 4 (1 + a:) - * dx = C K , m a n 



Thus, 15.11f) follows from the last two estimates. Using (|5.11|) it follows that 
>i>ast,i>a>B>t>) SHf(a', s',t') d(j,(a', s',t') 

Q12 



< C m a 



2JB(t ,ri) JO 
m-2 
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for all m > 4. Using the other cases of Lemma f5. 71 one can show similar estimates 
for the integrals over the sets Qio and Qn. This shows that SH<f,f 1 (a, s, t) decays 
rapidly for a — > uniformly over B(to,rj) x [—2, 2]. 
Let us consider now the term corresponding to fi- 

SH<pf 2 {a,s,t) = ((j>f2,ipast) = / (<t>i>ast,i>a's't')SHf(a',s',t')dLi(a',s',t'). 

JQi 

We will decompose Q 2 = S21UQ22, corresponding to ||(t — t')|| > rj and ||(t — t')\\ — 
77, respectively. Observe that, for \\(t — t')|| > rj and K > 1, 

JB{t ^Y \ a\ ) Jr, V a iJ \ a iJ 

Observe that, on the region Q 2 i, the function SHf(a', s', if) is bounded by C (a') 3 / 4 
since / is bounded. Thus 



/ {4>i>ast,^a> s <t>) SHf(a', s', t') dfj.(a', s' , t') 
JQ21 




and this is of rapid decay, as a — > 0, uniformly over Q2i- As for the region Q22, if 
t € B(to, rj) and || (t — t')\\ > 77, then t' G B(to,6) and thus the function SHf decays 
rapidly, for a — > 0, over this region. Repeating the analysis as in the case Q12, 
we can prove that jQ 22 (<fiipast,i>a's't') ^ f( a ' ',s' ,t')dfx(a' ,s' ,t') is of rapid decay, 
as a — > 0, uniformly over Q 2 2- Combining these observations, we conclude that 
SH<pf 2 (a, s, t) decays rapidly as a — > uniformly over £?(io, 7?) x [—2, 2]. 

It follows that SH g (a, s,t) decays rapidly as a ~> uniformly for all (s,t) E 
B(t , 77) x [—2, 2] and, thus, by Lemma l53l gi(^) decays rapid as |£| — * 00. We can 
now conclude that g decays rapidly as |£| —> 00, hence completing the proof of (i). 

For part (ii), we only sketch the idea of the proof, since it is very similar to part 
(i). Let (to, So) € T>. We consider separately the case |so| < 1 and |so| > 1. In the 
first case, for all t € B(t , 5) and s G B(s , A), we have that \SHf(a, s, t)\ = 0(a k ), 
as a — > 0, for all k € N with the 0(-)-term uniform over (t, s) e B(to, S) x B(sq, A). 
Choose 4> € L 2 (M. 2 ) which is supported in a ball B(tQ,v) with v <§; 6 and let 
77 = |. Then the proof proceeds as in part (i), replacing B(to,5) x [—2,2] with 
BitojS) x B(sq.A). Also, for the estimates involving inner products of ipast and 
ipa's't' we will now use Lemma |5.7I including the directionally sensitive term. For 
example, when i/a, \fa' < 5, by Lemma 15.71 we will use the estimate 

m.,„M\ < c K (1 + a)-" (1 + (1 + , 

rather than (j5~TTj|) . We can proceed similarly for the other estimates. The proof 
for the case |so| > 1 is exactly the same, with the transform SHf {a, s, t) replacing 
SH f (a,s,t). □ 
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6. Extensions and generalizations of the continuous shearlet 

transform 

As mentioned above, there are several variants and generalizations of the con- 
tinuous shear lets introduced in Section [3J In fact, using the theory of the affinc 
systems, we can obtain several other examples of continuous wavelets depending on 
the three variables: scale, shear and location. 

For example, we can generalize our construction by considering the case where 
A is given by i|3.1f) and M is of the form 

(a -a s s\ 
M s = = B As, aeI,seS, 

\0 a* ) 

where B = ( 1 s \ Ax — ( 1 and < S < 1 is fixed. If S — h, then we have 

Vow \o a v __ 2 

the case of continuous shcarlcts. In general, for other choices of 8, As will provide 

different kinds of anisotropic scaling. Using a construction similar to the continuous 

shearlets in Section |3 for each < 5 < 1, one can construct systems of the form 

{^ast=T t D Ms iP: (M 5 ,t) € A}, 

where tj) as t is supported on the set: 

supp^ ast C {(6,6) : 6 € U [i |], | fa - s\ < a 1 - 6 }. 

This provides a new family of transforms SH S j(a, s, t) = (f, ip a st), for various values 
of S. It turns out that, provided < 8 < 1, the behavior of the transforms 
is very similar to the continuous shearlet transform in dealing with 
pointwise and linear singularities. More precisely, one can repeat the analysis in 
Sections 14.11 14.21 and 14.31 using SH S j:(a, s,t) (for 6^1) rather than the continuous 
shearlet transform. Similarly, the behavior of the transforms SH S f(a, s,t) is very 
similar to the continuous shearlet transform in dealing with curvilinear singularities 
(see Section |4~4"|) . provided < S < 1. However, the proof of Theorem 15. II requires 
5 = h; that is, we need to use the continuous shearlet transform. 

Another variant of affine systems generated by A, given by i|3.1fl . is obtained by 
reversing the order of the shear and dilation matrices, namely, by letting 




= As B, ae I,s e S, 



where As, B are defined as above, and < 6 < 1 is fixed. Also in this case, we can 
construct variants of the continuous shearlets. Using the same ideas as above, we 
obtain a system of functions ip as t with support 

supp^C {(6,6) :6e [-!,-£] u[£,j[], l^^-sl^a 1 ' 5 }. 

It turns out (as one can see from the support condition) that the transform associ- 
ated with these systems is not even able to 'locate' the linear singularities, in the 
sense described in Section [Ol 

As we mentioned above, because the frequency support of the continuous shear- 
lets is symmetric with respect to the origin, the continuous shearlet transform is 
unable to distinguish the orientation associated with the angle 9 from the angle 
9 + 7r. In order to be able to distinguish these two directions, we can modify our 



28 



GITTA KUTYNIOK AND DEMETRIO LABATE 



construction as follows. Let ip € L 2 (R 2 ) be defined as in Section except that 
supp?/>i C [|,2]. That is, we have a one-sided version of the shear let ip illustrated 
in Figure 1. It is then clear that, if we consider the affine system generated by ip 
under the action of A given by ((3.1(1 and 1(3 .2(1 . this cannot provide a reproducing 
system for all of L 2 (R 2 ). In fact, ip is a wavelet for the subspace L 2 (H) V , where 
H = {(6,6) G M. 2 : 6 > 0}. In order to obtain a wavelet for the space L 2 (M 2 ), 
we can modify the set A as follows. Let A' given by 1)3.1(1 . where G c GL2(M.) is 
the set of matrices: 

( (ia -tjas\ 
(6.1) g=\m= , a e I, s e S, £ = -1, 1 

where / C K + , S C 1 Then the function ip is a continuous wavelet for L 2 (R 2 ) 
with respect to A'. These modified versions of the continuous shear lets depend on 
four variables a, s, t, £ and have frequency support: 



SUpp Ipastl C 



'{(6,6) : 6 e [-1 , -U || - *| < if ^ = -1, 



J(6,6):6G [i,f], ||-s|<^}, if* = i. 

We remark that these modified versions of the continuous shearlets are in fact 
complex functions, whereas the shearlets we use throughout this paper are real 
functions. 

Finally, there exists a natural way to construct continuous shearlets also in di- 
mensions larger than 2. We refer to ^5] for a discussion about the generalizations 
of shear matrices to higher dimensions. 

Appendix A. Additional Computations 
Proof of Theorem 12.11 

Suppose that 1)2.4(1 holds. Then, by applying Parseval and Plancherel formulas, 
for any / 6 L 2 (R n ) we have: 

\(f,T t D M ip)\ 2 dX(M)dt 



2 



det M\d\{M)dt 
\AetM\dX{M)dt 

2 

dt | det M\dX(M) 



(MM*-)) (t) 

(/V)(M*-)) V (t) 
l/(6| 2 |^(M t 6| 2 |detM|d^A(M) 



G JR" 



= / |/(OI 2 A(V0(£)^ = ||/|| 2 . 



Equation ((2.3(1 follows from the above equality by polarization. 
Conversely, suppose that 



\(f,T t D M iP)\ 2 dX(M)dt = \\f\\ 



R" JG 
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for all / € L 2 (R n ). Let £ be a point of differentiability of A(V0(£) and let /(£) = 
\B(£o, r)\~ x / 2 XB(£ ,r){Cji where B(^ ,r) is a ball centered at £ of radius r. By 
reversing the chain of equalities above we conclude that 

l-B(fo,r)| Js(c ,r) 

for all r > 0. Letting r — > 0, we obtain that A(-0)(£o) = 1, an d, since a.e. £ e M" 
is a point of differentiability, (|2 .4|> holds. □ 

The proof easily extends to function / S L 2 (U) V . In fact, it suffices to replace 
/ with fxv i n the argument above. 

Proof of equality (|3?7|) . 

Using Plancherel and Parscval formulas we have that 



\(f,T t W)\ 2 dt 



2 



dt 



(/w) (*) 



dt 



(A.l) 

Using a similar computation we have 

f'2 



i/(ohw(oi 2 

/ / I \{PcJ^ as t)f ^dsdt 
Jr J -2 Jo a 



-2 JO 
1 



(A.2) 



-2 JO 
2 



/(0xc 1 (C)V'(M* s C)e 2 "« t ^ 

2 



da 
7,3/2 



ds 



/ XCl ^(M* s -)) (*) 



, da 
at — ds 



2 Jo 



\fm 2 xcAO\4>(MU)\ 2 dt; 



1 3 / 2 

da 
375 



As in the proof of Proposition 13. 21 for £ 6 Ci, we have that 



'l^i^i)! 2 — - 



Thus, using the last equality, from (| A. 2|> we have that 



(A.3) 



da 



-2 Jo 

Similarly, we have that 



\(PcJ^ as t)V ~^dsdt= / 1/(01^(0 / \M<*i)f — dt 



da 



(A.4) 



2 Jo 



2 da 



■ ds dt 



IJWxc 2 (0 / hM< 2 )| 2 -d£ 



Thus, combining ()A.ljl . (|A.3fl and 1A.4|I and using (|3.6[) we have that 



\(f,T t W}\ 2 dt- 



I I I \{PcJA as t)\ 2 ^dsdt 

JR J — 2 J a 



c-2 
-2 JO 



(Pcj^ii) 



2 da 



ds dt 
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